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States supported by chaotic open quantum systems fall into two categories: a majority showing 
instantaneous ballistic decay, and a set of quantum resonances of classically vanishing support 
in phase space. We present a theory describing these structures within a unified semiclassical 
framework. Emphasis is put on the quantum diffraction mechanism which introduces an element of 
probability and is crucial for the formation of resonances. Our main result are boundary conditions 
on the semiclassical propagation along system trajectories. Depending on whether the trajectory 
propagation time is shorter or longer than the Ehrenfest time, these conditions describe deterministic 
escape, or probabilistic quantum decay. 
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Introduction. — Quantum states populating 'open' 
chaotic cavities decay to the outside environment and, 
thence, have the status of resonances. In spite of the 
ubiquity of the general setup — open quantum chaos 
is realized in many of the devices currently explored 
in mesoscopic physics, quantum optics, and cold atom 
physics — salient features of these resonances are not 
fully understood. While the deep quantum regime (the 
Ehrenfest time, ig, marking the diffractive disintegration 
of minimal wave packages shorter than classical escape 
times, td) appears to be under reasonable control [T], it 
is the opposite, semiclassical limit which poses unsettled 
issues [2]- 

Broadly speaking, the states populating an open cav- 
ity can be grouped into two families: states evolving near 
classically and escaping deterministically after a classical 
flight time, and a fraction ~ exp(— fg/td) of quantum 
resonances, whose probabilistic decay is characterized by 
a finite imaginary offset iT/2 to the real resonance energy 
E. The most basic quantity characterizing the statistics 
of resonances of complex energy z = E + i-^ is the res- 
onance density p(z). Although the quantitative profile 
of that quantity is not fully understood, the density ap- 
pears to be gapped against the real axis, T = (the 
existence of rare midgap states notwithstanding [21 [3].) 
The integrated number of resonances at a given value of 
E has been found to obey the so-called fractal Weyl law, 
p oc hr d s , where df is a non-universal fractal exponent. 

Previous work on the phenomenon includes the for- 
mulation of lower bounds on the resonance gap [IHS], 
semiclassical approaches based on short periodic orbits 
trapped in the open system [7j, a description in terms of 
non-unitarily evolving Husimi functions [5] , phenomenol- 
ogy based on a mixture of phase space dynamics and 
random matrix theory, resp. [Uj, and numerical analy- 
ses [21 [5l I8HT2]. However, a unified theory of resonance 
formation in terms of first principle semiclassical dynam- 
ics appears to be missing and the formulation of such a 



theory is the subject of the present work. 

Specifically, we will explore the quantum dynamics of 
states concentrated on classical trajectories in terms of 
phase space Wigner functions. Assuming globally hyper- 
bolic classical dynamics we will describe how quantum 
fluctuations operational on long trajectories convert the 
deterministic classical escape of short trajectories into 
probabilistic quantum decay. 

Our analysis is organized in three conceptual steps. 
We first consider a one-dimensional toy model to intro- 
duce the phase space language used in the rest of the pa- 
per. We then generalize to the more complex setting of a 
higher dimensional cavity, and derive effective boundary 
conditions determining the decay rates of the system. In 
a third step, we analyze these equations for both short 
and long trajectories. 

Toy model. — Consider a one-dimensional 'cavity' pa- 
rameterized by the spatial coordinate q € [— qo, qo], while 
coordinates to the right/left of +qo/—qo define connect- 
ing 'leads'. We assume free intra-cavity particle dynam- 
ics, Hq — p 2 /2m and, crucially, no backscattering barri- 
ers at the cavity/lead interfaces. 

Life times and energies of the resonant states sup- 
ported by the system may be calculated by matching 
solutions of the cavity Schrddinger equation to outgoing 
boundary conditions [T3J [T3] , i.e. by requiring that cavity 
wave functions tp(q) and their derivatives smoothly con- 
nect to vacuum wave functions ip± (q) = a± exp(±«<7/-c) 
at the right/left interface. Here, a± — const., and 
k = k — i^pisa complex wave vector whose real- and 
imaginary-part define the energy hk = (2mE) 1 ^ 2 and 
life time hk^ = T/v of resonant states, resp., where 
v = hk/m. The divergence of the reference states at 
spatial infinity q — > ±oo, is a formal means [13] to the 
fixation of decay rates, as exemplified below. 

For the intra-cavity wave function we make an ansatz 
V>(?) = E CT =± ae <T( ^ (9)+fcrt " (<?)/2) in terms of left- and 
right-propagating partial amplitudes where a = const., 
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and the somewhat unconventional denotation <fi(q) = kq 
and t(q) = q/v, for the real and imaginary contribution 
to the phase, resp., will be motivated shortly. With this 
choice, the boundary conditions obtained by matching 
wave functions and their derivatives at the left and right 
interface reduce to the single algebraic equation, 



-(i2<P(q )+k r vt(q )) _ ^9 



k-i\ (vd g t-l) 



0„ 
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(1) 



Before evaluating this equation, let us translate from 
the language of wave functions to a phase-space for- 
mulation. To this end, we introduce the Wigner func- 
tion W(q,p) = J(da) e'^ipiq - ^)ip{q + §), where 
(da) = da/(2irh). For our specific system, 



W{q,p) = a 2 8{p- ( jp{E))e' 7kTVt ^ 



(2) 



where, a = ± labels the Wigner transform of the left and 
right moving components, resp., and the ellipses denote 
rapidly oscillating interference contributions. In discard- 
ing the latter, we loose track of the global phase of the 
wave function, while the information on amplitudes and 
phase derivatives necessary to evaluate boundary condi- 
tions is retained. Indeed, it is straightforward to check 
that 



(3) 



a 2 e «krvt( q) = |^ (?) |2 = j {dp)Wa{q:p)i 
J(dp)pW <T (q,p) 
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For the simple Id system, the linear dependence <p q 
implies d q 4> — k so that reduces to 



8E 
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kq 



(4) 



where we have 'fixed a gauge' e~ 2l ^ qo ' ) = i for the arbi- 
trary phase of the wave function and neglected contribu- 
tions fcr/fc <C 1. 

For the toy model at hand, vd q t(q) = 1, which means 
that the right hand side of Q vanishes, and r — > oo is the 
only consistent solution. This reflects the fact that a wave 
function will 'decay' with probability unity upon passing 
the reflectionless boundaries of the system. We next dis- 
cuss how the situation changes upon generalization to a 
higher dimensional system with chaotic dynamics. 

Chaotic cavity. — We consider a two-dimensional cav- 
ity with ballistic Hamiltonian H = p 2 /2m and chaotic 
boundary scattering. The cavity is open such that after 
an average time td, much shorter than any of the relevant 
quantum time scales, trajectories escape through one or 
several reflectionless openings. We define the Wigner 
function of the system's resonance states by obvious gen- 
eralization of Eq. ([2]), i.e. W(c[,p) — J(d 2 a) e~ lE K~Tp(c[ — 



f )'0(q + f )• To obtain the intra cavity evolution equa- 
tions of W, one adds and subtracts the Schrodinger equa- 
tions of the resonances -0 and ip to obtain [15] 

[H * W] + =2EW, [H*W]_ = -iT W. (5) 

Here, H = p 2 /2m is the Hamilton function and 
[A * B] T = A * B =F B * A where the Moyal product 
of phase space functions A = A(q, p) is given by [16] 
A* B = AB + f{A, B} + 0{h 2 ), and { , } is the Poisson 
bracket. 

We next consider the vicinity of an exceptionally long 
trajectory 70 spending time T ^> td inside the cavity. 
For completeness we note that long trajectories in open 
systems are found with low probability ~ exp(— T/tj). 
They typically form in the phase space neighborhood of 
strange repellers realized through periodic orbits trapped 
in the interior of the cavity (see below for further com- 
ments on this aspect). Assuming global hyperbolicity 
of the dynamics, we introduce a trajectory coordinate, 
q G [— vT/2, vT/2], a conjugate momentum p — p(H) — 
(2mH) 1 / 2 transverse to the shell of conserved energy, and 
a pair it, s of locally unstable and stable coordinates. In 

the asymptotic neighborhood of 70, the Hamiltonian can 
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then be approximated as H ~ H = + Xus, where A 
is a Lyapunov exponent. The corresponding dynamics is 
generated by [H *,.]- = ih{H Q , } = —ih£, where the 
Liouvillian 



C = vd q + X(ud u — sd s ) 



(6) 



describes propagation in the direction of 70, and ex- 
ponential expansion/contraction in the u/s coordinate. 
Nonlinear corrections to Ho can be described as H = 
H n + V, where V — V{u, s) is a polynomial of degree > 2 
in the variables u, s. The corresponding modification of 
the dynamics, [V * ■]- = —ih(A£ + Q), comprises a 
weak alteration of the classical Liouvillian, AC, and a 
quantum generator 



Q= > r 

n+m> 1 



/ ^ n rn ''- ^ u s ? 



(7) 



where c nm = c nm (q,u, s) are coefficient functions whose 
detailed profile will not be of much importance through- 
out. Although both contributions are nominally small 
in u, s, the quantum generator Q, will be seen to have 
a regularizing effect on classical singularities [T7], which 
will ultimately shape the profile of the resonance density. 

Life time in a chaotic cavity. — Close to the trajec- 
tory, the first of Eqs. {5), [H * W} + ~ [H * W] + ~ 
2 (p 2 /2m) W — 2EW simply describes the on-shell fixa- 
tion p ~ (2mi?) 1 / 2 . Turning to the second equation, 



h (£ + AC + Q) W{q, u, s) = TW(q, u, s), 



(8) 



we first discuss the linear approximation, AC, Q = 0, 
before including the correction terms in a second step. 
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For AC, Q — 0, ([8]) becomes a first order differential 
equation which is solved in terms of a left- and a right- 
moving contribution 



W(q,p,u,s) = J2 a 2 8{ P ~ap{E))e 



(rkrvt(q,u,s) 



(9) 



cr=± 



structurally similar to Eq. |2|. Here, t(q,u,s) are ef- 
fective parameter functions generalizing t(q) of the toy 
model and evolving uniformly along the trajectories 7 = 
7 X piercing the phase-space point x = (q, u, s), Ct(x) = 1. 
To solve this (partial first order differential) equation, we 
consider its characteristics, i.e. the trajectory 7 X . On 
7 X , the equation assumes the form d r t(x(r)) = 1, where 
q( T ) — Q + v ft, u(t) = ue Xr , s(r) = se~ Ar and (q, u, s) 
are starting values of the evolution. We solve the char- 
acteristic equation as t(r) = t + t°, where r increases 
uniformly until 7 x (r) hits the effective boundaries of the 
problem, and t° is a freely adjustable parameter. 

To understand the role of the boundaries, note that 
7 X will leave the cavity either through a physical inter- 
face, S along with 70 = 7( g ,o,o) (cf. Fig. flj), or it will 
depart from 70 up to some classical threshold u ~ c 
within the cavity (solid line). We assume that points 
separated from 70 by scales ~ c have become generic 
and will exit in the classical, and hence negligibly short 
time ~ td- The union I = S a U C of S and the sur- 
face C = {u = c, s, q} then defines the effective 'out- 
going interface' of our problem. Similarly, the union 
lj = Si U Ci of the left vacuum interface Si and the 
surface C; = {u, s — c,q} defines the incoming inter- 
face. The travelling phase space point x(r) = (q,u, s)(t) 
hits the exit interface I , at the smaller of two times, 
t = t Q = t (q,u) = min(f - |,iln(j^)), depending 
on whether S or C Q is the terminal. Likewise, 7 X has 
entered the cavity through Ii at a large negative time 
t = -ti = -U{q,s) = -min(f + f , ^ ln( jfj )). Fixing the 
free parameter t° such that the temporal range of the tra- 
jectory is symmetric around zero, t(r = ti) = — t{r = t ), 
we find that the solution to Eq. ^ is governed by the 
function i(x) = \{ti{q, s) — t Q {q, u)) while T 7 = ti + t D is 
the intra cavity flight time of 7. Notice that for points 
x G I a at the exit interface, i (x) = 0, meaning that 
i(x) = T 7 /2 attains its maximal value. 

Finally, the boundary conditions Eq. are gener- 
alized by replacing the one-dimensional variable t(q) by 
t(q,u,s), and the derivative vd q by C, i.e. a derivative 
acting in the direction of the Hamiltonian flow [18] . The 
generalization of Eq. ([T]) thence reads 



8£; 



-rt/n _ 



= 1 - Ct, 



(10) 



where t — t(x) = T 7 /2, and x 6 7 is on the exit interface. 
Eq. ( 10 ) is a principal result of the present paper. In the 




FIG. 1: A phase-space point (q,u,s) in the vicinity of 70 
propagates along a unique classical trajectory, 7. It will exit 
the cavity either through the interface S or within the cavity 
through the surface C . Similarly, the union of the left vac- 
uum interface Si together with the manifold d defines the 
incoming interface. 



Short trajectories. — To start with, we consider trajec- 
tories 7 which are exceptionally long-lived, T = T 7 3> td, 
yet short in comparison to the scale tE ~ j ln(c 2 /H) 
where quantum uncertainty leads to the disintegration 
of semiclassically propagating wave packets. As long as 
T < tE, the quantum generator Q does not modify the 
dynamics in essential ways — a statement to be verified 
below — and the same goes for the generator of weak 
classical corrections, AC. We may thus take the bound- 



following we discuss its implications for different types of 
trajectories. 



ary condition ( 10 1 at face value, and conclude that due 
to the homogeneous evolution Ct = 1, the right hand side 
of the equation vanishes. As with the toy model, this im- 
plies a diverging decay rate, Y — > 00. In this divergence 
reflects the fact that wave packages travelling on classi- 
cally short trajectories leave the cavity with certainty at 
the exit point. 

Long trajectories. — We now turn to the case of long 
trajectories, T > ts- For asymptotically long trajec- 
tories, T — > 00, the classical shrinkage s — s(t) ~ 
cexp(— Xt) would lead to singularities in the function 
ti(q,s) and, ultimately, in the Wigner function. Within 
the present formalism, these singularities are regularized 
on time scales larger than ts, where s < s(ie) ~ h has 
shrunk down to quantum scales. On these scales, the 
quantum generator Q of Eq. Q is no longer small in 
comparison to the classical generator C (while the cor- 
rection AC to the classical flow continues to be largely 
irrelevant.) The ensuing modifications of the dynamics 
can be described in various ways (cf. Ref. 20 for a 
treatment tailored to the formalism applied here), the 
invariable conclusion being that the shrinkage of clas- 
sically evolving variables gets cut off by quantum fluc- 
tuations. Technically, this conclusion rests on the ob- 
servation that in the evolution equation for the vari- 
able s the higher order derivatives ~ 9™ >1 present in 
the quantum generator build up 'pressure' counter- 
acting the classical contraction. This is seen in explicit 
terms in the Fourier/Laplace representation of the evo- 
lution equation, where these derivatives assume the form 
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of algebraic factors, cutting the logarithmic 'ultraviolet' 
singularities of the classical equation. Referring to the 
supplementary material for more details [21], we note 
that to leading semiclassical accuracy functions which in 
the classical theory evolve as f(q, \u\, \s\) get replaced by 
f(q, \u\ + h/c, \s\ + h/c). Here, c is symbolic notation for 
classical (/i-independent) functions over which we have 
no explicit control, and the substitution \u\ — > \u\+h be- 
comes effectual in the large negative time asymptotics of 
a trajectory, where u rather than s scales to small values. 

To understand the consequences of this reguiarization 
mechanism, consider the trajectory time parameter, t = 
ti/2 at the exit point of 7. Now notice that ti(q, \s\) —> 
min(T/2 + q/v, A" 1 ln(|s| + h/c)) = A" 1 ln(|s| + h/c)) ~ 
tE, where we used that, T > ts- The crucial observa- 
tion here is that the reguiarization effectively truncates 
the in-time function t% at values tE- As a consequence, 
the interface derivative Ct = \C{ti — t Q ) = 1/2 reduces 
to one half of the value before quantum reguiarization. 
Substitution of this value into Eq. ( 10 1 shows that the 



quantum theory admits finite values of the decay con- 
stant, determined by 



ro 
2 



t E \ h 



h ( [8Et E 
V E ( l0g { — 



(11) 



where W is the Lambert function and ellipses denote sub- 
leading double-'log' contributions. Eq. (11) states the 



decay rate in terms of the Ehrenfest time in combination 
with non-universal short time cutoff h/E. However, in 
the semiclassical limit, h — > 0, the dependence on E drops 
out, and we are left with the asymptote Tq ~ hX. Before 
commenting on this result, we note that the appearance 
of a finite decay rate within our present formalism follows 
from the fact that, by Heisenberg uncertainty, quantum 
mechanics is not capable of resolving the phase space fine 
structures pertaining to the evolution of long trajectories 
T 7 > tE- Each such trajectory should, rather, be thought 
of as a distribution defined by the union of trajectories 
with uncertainty ~ h in their phase space coordinates. 
At a given instance of time, a fraction of this distribu- 
tion escapes, as described by the rate IV 

Effective decay rate. — Our above analysis was over- 
simplifying in that it treated escape from an isolated 
long trajectory, 70, as tantamount to escape into the lead 
vacuum. This picture ignores the fact that escaping tra- 
jectories may get 'folded back' into the repeller domain 
supporting 70, and thence be trapped again. The total 
escape rate thus presents itself as the rate of initial es- 
cape, ~ A -1 , [U[22], multiplied by the fraction of trajec- 
tories managing to escape the repeller, (1 — d), where d is 
the fractal (information) dimension d [H|24] of the latter. 
The effective decay rate, r — s- T = hX(l — d) is generally 
identified with the inverse of the classical escape time of 
the system. 



Fractal Weyl law and random matrix regime. — For 
completeness we note that a finite quantum mechanical 
decay rate is attributed to states located in the vicinity 
of exceptionally long trajectory. Due to the exponen- 
tial scarcity of these trajectories, the corresponding phase 
space measure scales as [TU] nr = £lEe~ tE ^ td ~ fi^/i 1 /*^, 
where £Ie is the phase space volume of the energy shell 
(in units h), implying that nr ~ hr di with fractal di- 
mension df = 1 — 1/Xtd [25]. We finally note that in 
the quantum regime, ig < td, which is complementary 
to the semiclassical regime td < tE studied here, random 
matrix scattering theory predicts [25] T ~ j^ln(Etd/h), 

where E is some cut-off energy scale. Comparison with 



(11) shows that the two results match at the boundary 
td ~ ts- However, we cannot say whether this matching 
is coincidental or not. 

Summary. — We have formulated a semiclassical the- 
ory quantitatively explaining the formation of quantum 
escape processes in open chaotic quantum system. Our 
theory describes the deterministic nature of escape on 
short trajectories, in contrast to the formation of finite 
decay rates by quantum fluctuations on long trajectories. 
These fluctuations effectively truncate the buildup of 
classical singularities in the system's hyperbolic dynam- 
ics. By the same token, they render the escape dynam- 



ics effectively stochastic, as described by the rate (11) 



Somewhat counterintuitively, the relevant escape dynam- 
ics takes place on long trajectories T > tE, while the 
ensuing decay rates are classically short V ~ h\(l — d). 
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We here discuss how quantum fluctuations regularize 
the unlimited classical contraction of the stable coordi- 
nate s in a system with globally hyperbolic dynamics. In 
the language of Eqs. (6) and (7) of the main text, the 
dynamics of the variable s is described by a differential 
equation of the structure 

sd s + £ c n h 2n+1 d 2 s n+1 I f(a) = -af(s) (S12) 

V n>l / 

where a > 0, and in a manner inessential to the present 
argument the coefficients c„ may depend on the variables 
q,u, s. 

Considering positive starting values, s > (the exten- 
sion to negative values is straightforward), we introduce 
a Laplace representation 



setting all c n = 0. The resulting function 



co 



dz i 



Co 



(S16) 



then displays the singular at small values of s plaguing 
the classical evolution equation of the stable coordinate. 



In the full solution Eq. (S15) the exponential factor 
e~^-"Si 2n+i ( riz ) cu ts the small s/large z singularity 
at values z ~ 1/h. The resulting integral can be esti- 
mated by a regularized function 



/(«) = 



co 



(s + ny 



(S17) 



f(s)= / dze- sz g{z) 
Jo 



(S13) 



in which (S12) takes the form pQ 



d z g{z) = - 



1 -a 



Y^c n h 2n+1 z 2n \g(z), (S14) 



The general solution of this equation is found by straight- 



forward integration over z, and when inserted into (S13) 
gives 

/•OO 

f(s) = c dze-^z^e-^i^^-) 2 ^ 1 ( S15 ) 
Jo 



with an integration constant cq. Eq. ( S15 ) now illustrates 



the role played by higher differential operators in (S12). 



To make the point, let us for the moment consider the 



first order differential equation obtained from (S12| by 



Finally notice that our argument crucially relies on as- 
sumed positivity of the coefficients c„. While the present 
construction cannot prove this feature, positivity is re- 
quired on principal grounds to ensure stability of the dy- 
namics. (Otherwise the Wigner distribution would cease 
to exist.) To actually demonstrate this stability, one has 
to work harder as in, e.g., Refs. A discussion tai- 

lored to the present formalism is contained in Ref. [4]. 
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